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ABSTRACT 

Aims. The DQ white dwarfs are stars whose atmosphere is enriched with carbon, which for cool stars (r eff < 8000 K) is indicated 
by the Swan bands of C2 in the optical part of their spectra. With decreasing effective temperature these molecular bands undergo a 
significant blueshift (~ 100 - 300A). The origin of this phenomenon has been disputed over the last two decades and has remained 
unknown. We attempt to address this problem by investigating the impact of dense helium on the spectroscopic properties of molecular 
carbon, the electronic Swan band transition energy T c and the vibrational frequency w e , under the physical conditions encountered 
inside helium-rich, fluid-like atmospheres of cool DQ white dwarfs. 

Methods. In our investigation we use a density functional theory based quantum mechanical approach. 

Results. The electronic transition energy T e increases monotonically with the helium density (Ar c (eV) ~ 1.6 p (g/cm 3 )). This causes 
the Swan absorption to occur at shorter wavelengths compared with unperturbed C2. On the other hand the pressure-induced increase 
in the vibrational frequency is insufficient to account for the observed Swan bands shifts. Our findings are in line with the shape 
of the distorted molecular bands observed in DQp stars, but the predicted photospheric density required to reproduce these spectral 
features is one order of magnitude lower than the one predicted by the current models. This indicates pollution by hydrogen or reflects 
incomplete knowledge of the properties of fluid-like atmospheres of these stars. 

Conclusions. Our work shows that at the physical conditions encountered in the fluid-like atmospheres of cool DQ white dwarfs the 
strong interactions between C2 and helium atoms cause an increase in T e , which should produce a blueward shift of the Swan bands. 
This is consistent with the observations and indicates that the observed Swan-like molecular bands are most likely the pressure-shifted 
bands of C2 

Key words, atomic processes - dense matter - stars: atmospheres - stars: white dwarfs 



1. Introduction 

A DQ white dwarf is a star with a helium-dominated atmosphere 
enriched with carbon, which is believed to be dr edged-up from 
the c ore by the deep helium convection zone (Pell etier et all 
1986). In the cool DQ white dwarfs <T eff < 8000 K) car- 
bon shows its presence by the C2 Swan bands in the optical 
spectrum. The two until now unexplained phenomena are ob- 
served in the local samples of these stars. There is a cutoff in 
their cooling sequence at r e ff ~ 6000 K (iDufour et all [2005; 
iKoester & Knist). |2008l) and an appea rance of "peculiar" DQ 
stars (following lHall & Maxwelll (12008) we call them DQp white 
dwarfs) at lower temperatures. These peculiar stars show Swan- 
like bands that are blueshifted by 100 to 300 A (Schmidt et al., 
1 1 9951; iBergeron et all \l99% lHall & Maxwelll 120081) . The repre- 
sentative spectra of normal and peculiar cool DQ stars are given 
in Fig . [D In the detailed analysis of D Qp stars ISchmidt et al.l 
dl995l 119991) and lHall & Maxwelll (120081) doubt the possibility 
that the shifts of the Swan bands observed in spectra of these 
stars are solely due to a strong magnetic field, which is not 
observed to be stronger than in "normal" DQ stars. Observing 
identi cal shifts of t he bands in all at that time known DQp 
stars, ISchmidt et al.l (1 19951) concluded that a different molecu- 
lar species could produce the observed bands. Analyzing the 
abundances of chemical species in the H/He/C mixture under 
the physical conditions found in the atmospheres of helium-rich 
stars, they proposed the C2H molecule as a source of the opti- 
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cal spectral features in the DQp stars. The idea that hydrogen is 
responsible for the spectra of DQp stars is consistent with two 
findings regarding the chemical evolution of cool white dwarf 
atmospheres. The already mentioned observation of DQ— >DQp 
transition indicates some sort of transformation of the physical 
or chemical properties of the atmospheres of these stars. Also the 
analyses of the local populations of white dwarfs by improved 
atmosphere models, which account for physics and chemistry 
of dense media, show the decrease in number of the helium- 
rich a tmosphere stars with T e g < 60 00 K (IBergeron et all 
I1997L 12001b iKowalski & Saumonl 120061) and the appearance 
of w hite dwarfs with atmospher e s highly enriched in hyd ro- 
gen (Kowalski & Saumo E 120061 iKilic et all H008, 2Q09aUbb. 
Subsequent work on the coolest stars indicates that the stars at 
the end of the white dwarf cooling sequence posses pu re hy- 
drogen atmospheres dHall et all 120081; iKilic et all l2Q10t) . Both 
these findings indicate that the accretion of hydrogen from the 
interstellar medium may change the composition of a helium- 
dominated atmosphere into a hydrogen-rich one, and shows that 
the appearance of DQp stars may be someh ow connected to 
this pr ocess. However, in the recent work of lHall & Maxweil 
d2008l) the authors ruled out the possibility that C2H molecule 
or other carbon and hydrogen bearing species are responsible 
for the absorption in the optical spectra of DQp stars and con- 
cluded that C2 itself is a reasonable candidate for producing the 
observed features. This claim is supported by matching of the 
long-wavelength edges of th e observed bands in DQ and DQp 
stars, as indicated in Fig. Q] Hall & Maxwe ll (2008) postulate 
that the Swan bands distortion could be produced by the high 
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Fig. 1. Optical spectra of the white dwar fs LHS179 (DQ, uppe r 
panel) and LHS290 (DQp, lower panel) dBergeron et alll 19971) . 
The solid line represents the fit to the spectrum of LHS179. 
The parameters of the fits are: T e s = 6500 K, log g=8 (cgs), 
and log C/He = -6.4. These values are in line with the atmo- 
spheric parameters derived in iDufour et al.l d2005l) . The vertical 
dotted lines indicate the positions of the Swan bandheads for 
Av = +2, +1, 0, - 1 from left to right respectively. The synthetic 
spectrum matches the entire known spectral ene rgy distribution 
of the star (given by BVRIJHK photometry of Ber geron et alj 
(Il997ln . 



pressure or the high excitation of the rotational states of C2. This 
claim is however purely speculative and is not based on strong 
theoretical or experimental evidence. 

The majority of DQ star atmospheres are highly depleted 
with hydrogen, and the dominant constituent species is he- 
lium. With cooling, an atmosphere of a helium white dwarf 
becomes more charge-neutral, less opaque and more dense at 
the photosphere. Also the amount of carbon, which signifi- 
cantly contributes to the opacity by deliv ering free el e ctrons , 
diminishes with the eff ective temperature dDufour et al.l (|2005); 
iKoester & Knisj d2008l) ). causing an increase in the photospheric 
density. The photospheric densities as a function of effective 
temperature for DQ white dwarfs with representative amounts of 
carbon are given in Fig. [2] The atmospheres of these cool stars 
(Jeff < 8000 K) posses extreme densities (as high as few g/cm 3 ) 
and represent a dense fluid, in which stron g inter-particle inter- 
actions affect the chemistry and physic s dBergeron et all f l995: 
Kowalski & S aumonl l2004t iKowalskil 120061: iKowalski et all 



20071) . Plotting the photospheric densities of a DQ star cooling 



sequence i n Fig. [H which is based on the carbon abundances 
derived by iDufour et al.l d2005l) . we see that the atmospheres of 
DQ stars reach fluid-like densities and therefore high pressures 
at the temperatures of the observed DQ— >DQp transition. This 
suggests the pressure effects as a good candidate to be responsi- 
ble for the spectral distortions observed in DQp stars. 

In order to understand the origin of the spectral features 
observed in DQp stars, we investigate the impact of dense he- 
lium fluid on the spectroscopically important properties of C2 
molecule. We performed ab initio calculations of C2 in dense 
helium, looking for changes in quantities that shape the optical 
spectrum of C2, the electronic transition energy T e , and the vi- 
brational frequency cj e . 



2. Computational approach 

The most widely used method for the quantum mechanical treat- 
ment of dense, extende d many-particle systems is the density 
functional theory (DFT, iKoch & Holthausenl d2000l) ). Over the 
last decade DFT methods proved to be very useful in calcu- 



lations of various chemical and physical properties of atoms, 
molecules, solid states and matter under high compression. In 
our work we used one of the most common DFT version, the 
generalized gradient app roximat i on (G GA) with PBE exchange- 
correlation functional (Perdewl 119961) . The calculations were 
perfo rmed by the plane-wave DFT ABINIT code ( Gonze et al., 
120021) with the ultra-soft pseudopotentials dVanderbiltt 11990) 
(PAW method and the energy cutoff of 15 a.u.). In order to ob- 
tain the distributions of helium atoms around C2 we derived the 
He - C2 interaction potential using a combined DFT, molecular 
dynamics, and classical theory of fl uids approach. First we p er- 
formed Car-Parrinello simulations dCar & Parrinellol 1 19851) of 
C2 emerged in dense helium at density of 0.5 g/cm 3 and the tem- 
perature of 5802 K (equiv. of 0.5 eV). The periodically repeated 
simulation box contained one C2 molecule and 64 He atoms. 
As the ABINIT code is not well suited for t hese calculations, 
they were performed with th e CPMD code dCar & Parrinellol 
1985| iMarx & Hutterl 120001) . keeping the same psuedopoten- 
tials, DFT functional and energy cutoff. We checked that us- 
ing the same computational setups both DFT codes give the 
same answers for the total energies (within 10~ 5 a.u.). A trajec- 
tory of one million configurations was generated (lOOps long). 
The C2 - He interaction potential was derived by solving the 
Ornstein-Zernike equation in the Percus-Yevick approximation 
dMartvnovlll992l) . taking as input the pair-distribution function 
obtained from the simulation. In order to calculate the energies 
of C2 in dense helium in 3 n„ and 3 H g states, which constitute 
the Swan band transition, for a given T and p we calculated 
the ground state energies of a set of 50 ionic configurations 
obtained from the classical Monte Carlo sim ulations of helium 
fluid w ith one C2 molecule, performed using Al len & Tild eslev 
(1987) Monte Carlo code with the der ived C2 - He i nterac - 
tion potential and He - He interaction of lRoss & Young| (fl986). 
We used this approach because an explicit quantum molec- 
ular dynamics (QMD) simulations would be computationally 
much more demanding, and be undoable for extended systems 
with hundreds of atoms. The relevant 3 n„ and 3 Tl g states are 
constructed according to the electronic configurations given by 
Balli k & Ramsavl (Q963). The symmetries of the five lowest KS 
orbitals of C2 molecule are: o- g (s ),a u {s),n u (_p),n u (p),o- g (p), in 
agreement with Ballik & Rams a v d!963h . The lower 3 n„ state 
((cr g 2s) 2 (cr u 2s) 2 (n u 2pj i (o- g 2p)) is obtained by fixing the occu- 
pation of spin-up orbitals as 11100 and spin-down 11111. The 
upper 3 n g state ((cr g 2s) 2 (cr u 2s)(7T u 2p) 3 (cr g 2p) 2 ) is obtained by 
fixing the occupation of spin-up orbitals as 10101 and spin- 
down 11111. Then the Swan Band electronic transition energy 
is given as the difference of these two calculations. The obtained 
equilibrium separations for these two states, 2.49 a.u., 2.37 a.u. 
and the electronic transition energy T e = 2.492 eV agree well 
with the ex perimental values ( 2.48 a .u., 2.39 a.u. and 2.39 eV re- 
spectively, iBallik & Ramsavl d 19631) ). which justifies the usage 
of DFT. 

The impact of dense helium on the vibrational frequencies 
of C2 was investigated by computing the vibrational density of 
states by the Fourier tran sform of the velocity autocorrelation 
function dEgelstaffl. [1991 derived from the velocities obtained 
from Car-Parrinello molecular dynamic runs for selected densi- 
ties and C2 in a ground 'X state. 

In order to compute the atmosphere models we used our own 
stellar atmosphere code that accounts for various dense medium 
effects such as refraction, the non-ideal equation of state, and 
chemical equilibriu m. The intern al partition functions of carbon 
molecule is that of llrwinl d!981l) . For the carbon atom and ion 
we set the T — OK values. This approximation is justified be- 
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Fig. 2. Left panel: Density at the photosphere (the Rosseland 
mean optical depth r r = 2/3, solid lines) and at the optical depth 
of the Av = Swan band formation region (t /1=5100 ^ = 2/3, dot- 
ted lines) in DQ white dwarfs of various He/C atmospheric com- 
positions: log He /C = -4, -5, -6, -7 and pure-He from bottom 
to top. The dashed line represents the photospheric density of the 
DQ star c ooling sequence, a ssuming the carbon abundance de- 
rived by iDufour et al.l (l2005h . No hydrogen bearing species are 
present in the models. Right panel: The change in the electronic 
transition energy AT e for the Swan bands transition in C2 as a 
function of helium density. Error bars represent the results of 
calculations fixing inter-atomic distance in C2 at values derived 
for isolated C2. The square at p — 0.25 g/cm 3 shows the result 
of the same calculation, but performed on the distributions of 
He atoms obtained from a quantum molecular dynamics simu- 
lation. The diamonds represent the results obtained relaxing the 
inter-atomic distance in C2 to the equilibrium position. The tem- 
perature is T = 5802 K = 0.5 eV. 

cause the electronic excitation energies of both carbon species 
are high, and for temperatures of a few thousand degrees the 
populations of the electronic excited levels are extremely small. 
Hydrocarbon species are not considered. 

3. Results and discussion 

It is well known that in helium fluid the interaction energies at 
the considered densities can be significant and may reach val- 
ues of ~ eVs (for instance the internal energy of atoms in dense 
helium is 0.5 eV/atom at 1 g/cm 3 and T = 5802 K). The in- 
teraction energy should in principle be higher for the species in 
the higher excited states, because their electronic charge is lo- 
cated farther away from the center of the atom/molecule. This 
should lead to an increase in an electronic transition energy. We 
computed the electronic transition energy of the Swan transi- 
tion in C2 that emerged in dense He at different densities and 
representative for the atmospheres of DQp stars temperature 
T = 5802 K = 0.5 eV. The obtained shift in the electronic transi- 
tion energy, AT e is given in Fig. [2] With the increase in the den- 
sity of helium, T e increases and the increase for densities up to 
0.25 g/cm 3 is linear and given by AT e (eV) = 1.6pH e (g/cm 3 ). 
This increase shows the smaller slope at densities between 0.5 
and 1 g/cm 3 and a similar increase at higher densities. The 
increase in the electronic transition energy should produce a 
blueshift in the Swan band spectrum, which is qualitatively 
consistent with the shifts observed for DQp stars. Positions of 
the band minima in the spectra of DQp stars imply a shift of 



Fig. 3. Left panel: The vibrational density of states of C2 at dif- 
ferent He densities. Right panel: the inter-particle distance in C2 
as a function of helium density (solid line) obtained for l X state 
and the average inter-particle separation in dense helium (dotted 
line). 



~ 0.08 eV. This indicates that the bands form at a density of 
~ 0.05 g/cm 3 , which is one order of magnitude smaller than 
the predicted photospheric densities in the models of these stars 
(~ 0.4 g/cm 3 , Fig. H). 

On the other hand the change in the vibrational frequency is 
prominent, but appears at significantly higher densities. The re- 
sult in terms of the vibrational density of states of C2 in 'X state 
is given in Fig. [3] The vibrational frequency is barely affected up 
to a density of 0.5 g/cm 3 and shifts toward higher frequencies at 
higher densities. This can be explained looking at the impact of 
compression on the mean inter-atomic separation in C2. In Fig. [3] 
we present the variation of this value as a function of density to- 
gether with the mean inter-particle separation in the fluid helium. 
The distance between carbon atoms that constitute the molecule 
remains unaffected up to a density of 0.5 g/cm 3 , and decreases 
slowly at higher densities. The decrease is correlated with the 
mean inter-particle distance in fluid helium. The weak density- 
dependence of the inter-particle distance in C2 (for considered 
densities) justifies the usage of fixed inter-particle distance in 
the C2 molecule in the calculations of T e . We also notice that 
to reproduce the observed shifts by the increase in the vibra- 
tional frequency, the difference in this increase should be at least 
1000 cm -1 (AE = 0.5^(A<y e U pp er - Aw e j ower )), which could oc- 
cur at densities hig her than 2 g/cm 3 . Potent ial rotational excita- 
tion, as proposed by Hall & Maxwelll (12008) for the reason of the 
shifts, is difficult to explain, because this is a temperature effect 
and should be stronger at higher temperatures, which is opposite 
to the temperature trend of the DQ— >DQp transition. Therefore, 
we conclude that the dynamical motion of carbon atoms is un- 
likely to be responsible for the distortion of Swan bands. Our 
results show that the distortion of the C2 bands in DQp stars is 
caused by the pressure-induced increase in T e . 

After computing the shifts in the electronic transition en- 
ergy for the Swan band transition we attempted to reproduce 
the spectrum of a DQp star. We did this by using the standard 
Swan band spectrum, shifted at a given atmospheric level by the 
derived value of Ar e . We note that although this is a quanti- 
tatively meaningful approximation, its validation would require 
the modeling of the absorption by C2 in dense helium, which is 
a complex task, restricted by the limited applicability of quan- 



4 



Kowalski: C2 in cool DQ white dwarfs 



T 






4000 4500 5000 5500 6000 

MA) 

Fig. 4. Fits to the optical spectrum of the DQp white dwarf 
LHS290 (thick solid line). All the fits are with T eS = 5800 K and 
log g=8 (cgs). The parameters of the fits are: H/He=0, C/He 
1 .25 • 10~ 7 , AT e = OeV and AT e = 1.6p He eV (solid lines, Mi 
and M 2 ); H/He=0, C/He = 1.25 ■ 10~ 7 , AT e = 0.2 p He eV (dotted 
line, M3); mixed composition He/H model, AT e = 1.6pH e eV, 
C/He = 10~ 6 and H/He = 6.75 • 10~ 3 (dashed line, M 4 ). All 
synthetic spectra well match the entire known spectral energy 
distribution of the star, as shown in Fig. [5] 

turn methods beyond DFT to many particle systems. In Fig. 
|4] we sho w the optical spectrum of the cool DQp white dwarf 
LHS290 (Ber geron et all 1 19971) together with a set of synthetic 
spectra. The overall spectral energy distribution of that star is 
best reproduced by models with T e ff = 5800 K (Fig. [5), and we 
assume this temperature in our analysis and fix the gravity at 
log g = 8 (cgs). The C/He and H/He abundances are fitted to best 
reproduce the peaks of the Av = 0, -1 bands. Assuming pure He 
atmosphere the strength of the Swan bands is reproduced with 
C/He = 1.25 ■ 10 , but with the computed correction for T e 
the spectrum is far too distortecQ. The observed spectrum can 
be fairly well reproduced by a model with photospheric density 
~ 0.05 g/cm 3 , or assuming that the T e dependence on the density 
is weaker, AT e ~ 0.2pe e - The first case is realized by the addition 
of hydrogen to the atmosphere, which increases the opacity and 
lowers the photospheric density. The required amount of hydro- 
gen is H/He = 6.75 • 10 3 . In both cases the observed spectrum is 
fairly well reproduced. The minima of the bands are blueshifted, 
as most of the absorption occurs close to the photosphere (Fig. 
[5). The long-wavelength parts of the bands resemble those of 
the standard Swan absorption, because part of the absorption oc- 
curs in the less dense upper atmospheric layers, where C2 is un- 
perturbed. We notice that the abundances of molecular carbon 
and the resulting strength of its molecular bands could also be 
affected by high density, which could eventually impact the re- 
ported carbon abundances. 

4. Conclusions 

The so-called DQp stars represent a puzzle in the understand- 
ing of evolution of cool, helium-dominated white dwarf atmo- 
spheres. The DQ stars disappear at T e s ~ 6000 K, and few stars 
with apparently distorted Swan bands were detected at lower ef- 
fective temperatures. All explanation through the formation of 
different species, like C2H, magnetic fields, or roto-vibrational 
excitations failed to explain the spectral features of these stars or 
definitely assign them as the distorted bands of C2. We show that 

1 We note that i t rese mbles the spectrum observed for a star G240-72 
in Bereero n et all dl997|). Fig. 30, which op ens the possibility that this 
star (Tj-jf = 5640 K. fBergeronTt al. (1997)) is not only magnetic, but 
also shows the extremely pressure-shifted Swan bands. 
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Fig. 5. Left panel: Fits to the entire spectral energy distribution 
of the DQp white dwarf LHS29 (thick solid line) given by 
BVIJHK photometry of [Bergeron etalj d2001l) (bars). The pa- 
rameters of the fits are the same as in Fig.[4j Right panel: Density 
profiles of the discussed models. The vertical dotted line marks 
the photosphere (r r = 2/3). Circles indicate the optical depth 
at which Av = band forms in models Mi, M3 (upper curve) 
and M4 (lower curve). The square indicates the optical depth at 
which distorted bands form in model M2. 

the distortion of Swan bands originates in the pressure-induced 
increase in the electronic transition energy between states in- 
volved in the transition. This results in a blueshift of the molec- 
ular bands minima, and explains why the red edges of the bands 
match the spectra of normal DQ stars (Swan bands). Our results, 
when applied to the current atmosphere models, predict Swan 
bands shifts that are too large compared with the observed ones. 
This indicates that the density at the photosphere of DQp stars 
does not excess 0.05 g/cm 3 , and the input physics in the models 
or the understanding of the atmospheres of these stars, especially 
the pollution by hydrogen, requires further improvements. 
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